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Abstract 

The complete charge-current density and field strength of an arbitrarily 
accelerated relativistic point-charge are explicitly calculated. The current 
density includes, apart from the well-established three-dimensional delta- 
function which is sufficient for its global conservation, additional delta- 
contributions depending on the second and third proper-time derivatives of 
the position, which are necessary for its local conservation as required by 
the internal consistency of classical electrodynamics which implies that local 
charge-conservation is an identity. Similarly, the field strength includes an 
additional delta-contribution which is necessary for obtaining this result. 
The Lienard-Wiechert field and charge-current density must therefore be 
interpreted as nonlinear generalized functions, i.e., not just as distributions, 
even though only linear operations are necessary to verify charge-current 
conservation. The four-potential from which this field and the conserved 
charge-current density derive is found to be unique in the sense that it is 
the only one reducing to an invariant scalar function in the instantaneous 
rest frame of the point-charge that leads to a point-like locally-conserved 
charge-current density. 



03.50.De Classical electromagnetism, Maxwell equations 



1 Introduction 

Local conservation of the charge-current density four-vector J M is a necessary 
consequence of the antisymmetry of the electromagnetic field strength tensor F^ v . 
Indeed, the identity 

3% = 0, (1.1) 
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derives from taking the four-divergence of both sides of Maxwell's inhomogeneous 
equation 

<9^ = -4ttJ m , (1.2) 

where the left-hand side vanishes after contraction because F^ v is antisymmet- 
ric. The vanishing of the four-divergence of the current-density is thus not an 
ordinary conservation law (i.e., a 'weak' conservation law subject to the field 
equations being satisfied) but an identity (i.e., a so-called strong conservation law) 
which has to be satisfied even if is a distribution rather than a smooth func- 
tion. This conclusion is absolutely general and should therefore be true for an 
arbitrarily moving relativistic point-charge, that is for the charge-current density 
of the Lienard-Wiechert field, which turns out not to be the case for the customary 
formulation of this field. 

In this paper we show that if the Lienard-Wiechert current is properly calcu- 
lated, which implies that the Lienard-Wiechert field strength must be supplemented 
by an additional 5-function-like field, local charge conservation is restored. This 
conclusion is obtained by using only well known physical concepts, and a few 
basic results of distribution theory, but at the expense of some lengthy calculations 
whose details could not be given in the letter: The locally-conserved current of 
the Lienard-Wiechert field, reference d). 

For ease of calculation we use the biquaternion formalism. While this choice 
is a matter of convenience for the present paper, it is a necessity for subsequent 
publications because the spinor decomposition of four- vectors enables complicated 
four-dimensional integrations to be made very efficiently and in full generality. 

The general structure of this paper is the same as the letter [1], in which the 
tensor formalism was used throughout. Thus we begin in section|2]by setting forth 
our notations and definitions, including a number of basic quaternion definitions 
to make our paper self-contained. A number of four-dimensional geometrical and 
kinematical identities needed in this and subsequent papers are derived, or else 
sufficient information is given to explain how they were derived. 

In section |3] we recall the customary formulation of the Lienard-Wiechert 
potential and field, and show that the corresponding charge-current density is not 
conserved. This requires calculations that are repeated in more details in the 
following sections, where it is also explained why the standard formulation does 
not lead to an identically-conserved current. 

The locally conserved charge-current density is derived in section HI The 
method used consists of postulating a very general form for the four-potential, 
expressed in the causal coordinate system defined in sectional and to seek under 
which conditions it leads to a conserved point-like charge-current density. This 
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current density, i.e., 

(1.3) 

T = T r 

has a rather complicated structure: Kinematically, it depends on the three invariants 
£, k, and as well as on the two four-vectors Z^ and Z^, and on the biacceleration 
Z^ through the invariant x\ geometrically, on the angular variables through the 
null four- vector K^(9, (f>); and, distributionally, on Dirac's 5-function and its first 
two derivatives. 

In section H] a straightforward derivation of the conserved charge-current den- 
sity is given. This mathematically rigorous derivation is based on a full exploita- 
tion of Schwartz's structure theorem of distribution theory, which explains how 
^-functions arise as partial derivatives of continuous functions [!2l |3). It will be 
seen that all is needed is a correct characterization of the nature of the singular- 
ity associated with a spin-less classical electron (namely a point-particle whose 
four-potential reduces to a scalar function in its instantaneous rest frame), and 
everything follows from a consistent application of Maxwell's theory. This char- 
acterization consists of replacing the Lienard-Wiechert potential by the expression 

(1.4) 

T=T r 

where the generalized function T(£) explicitly specifies how to consistently dif- 
ferentiate at £ = 0. Similar forms of the potential of a classical point-charge have 
been discovered independently by Frank R. Tangherlini [7], as well as by others 
flU, in various contexts. It then turns out that by starting from (11.41) to calculate 
the field F^ u and current J M according to their standard definitions one imme- 
diately finds the locally conserved current (11.31) . provided all 5-terms generated 
by differentiation are retained until the end of the calculation — that is provided 
Afj,, F^ u and J M are interpreted not just as distributions but as nonlinear generalized 
functions JH[6]|. 

In the conclusion, section [61 we discuss the consistency of our main results, 
equations (|1.3I) and (|1.4I) . with the standard formulations of Coulomb's potential, 
Lagrange's function, and Green's function of classical electrodynamics. Finally, 
we summarize these results in the form of a theorem, stating that the potential (11.41) 
is unique in the sense that it is the only four-potential reducing to an invariant scalar 
function in the instantaneous rest frame that leads to a point-like locally-conserved 
charge-current density. 



A^ = e-^T(0 
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2 Notations and definitions 



In this and related papers Hamilton's biquaternion (i.e., complexified quaternion) 
formalism is used instead of vector or tensor calculus. This is because biquater- 
nions provide concise and fully-general explicit expressions for all covariant quan- 
tities appearing in electrodynamics, including spinor and null quantities for which 
there are no explicit formula in ordinary vector or tensor calculus, e.g., equations 
(|2.1 II) and (|2.12l) below. This enables all calculations to be made explicitly and ex- 
actly (that is to obtain final results in closed form without using any approximation 
or limiting process) and within a fully-general and consistent four-dimensional 
framework that is eminently suited to the description of non-inertial motion in 
special-relativity. 

The same is true with regards to Cartan's differential forms, a formalism that 
is particularly well suited to abstract manipulations such as proving theorems, but 
not to explicit calculations. For instance, with respect to equation (11.11) . the general 
form of the divergence is V • J = *d * J where * is the Hodge dual — whereas, 
from a strict mathematical point of view, V • J = d^J^ only for a Cartesian basis 
(even though it is commonly understood by physicists that this expression has to 
be properly adapted when working in a curvilinear coordinate system as in the 
present paper). In response to this the biquaternion formalism provides, for the 
divergence in the four-dimensional space of special relativity, the general form 
used in equation (|2.3I) below, i.e., V o J, which in contradistinction to Cartan's 
expression *d * J is explicit and directly amenable to calculations. Indeed, in 
Cartesian coordinates Hamilton's operator V is given by equation (12.21) below, and 
in the curvilinear coordinate system used in the present paper by equation (|2.27l) . 

A primer of quaternion methods in physics is given in reference HI. However, 
as a quick start, we recall that the non-commutative algebra of quaternions consists 
of four-dimensional objects Q = s + v, were s and v are the scalar and respectively 
vector parts of Q. The quaternion product is defined as 

(a + a) (b + b) := ab — a ■ b + ab + ba + a x b, (2. 1) 

but detailed calculations requiring this formula are rarely needed since the crucial 
advantage of quaternions is that they can be manipulated as whole symbols, there- 
fore dispensing of complicated calculations at the component level, as in tensor 
and vector calculus. This is particularly useful in relativity and electrodynamics 
where the Minkowski metric is obtained by allowing quaternion components to be 
complex numbers. (Hamilton introduced the prefix 'bi' to denoted complexified 
objects, e.g., biquaternion, bivector, biscalar for complex number, etc.) 

For example, the four-gradient operator (i.e., the four-dimensional generaliza- 
tion of Hamilton's 'nabla') and the four-vector current-density in equation (11.11) 
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correspond to the biquaternions 

d 

V = — + V, and J = p-ij. (2.2) 

(We use Gaussian units and put the velocity of light in the vacuum equal to one.) 
Equation (11.11) then becomes 

Voj = 0, (2.3) 

where the unary operator ( ) denotes quaternion conjugation, namely the operation 
(s + v) = s — v such that QR = R Q, and the binary operator o as in Q o R 
means taking the scalar part of the product QR. As is shown in (8), products 
such as QRST..., where four-vectors and their quaternion conjugate alternate, are 
automatically covariant. For example, V o J is an invariant scalar. Moreover, for 
any four- vector Q, the product QQ = QQ is an invariant scalar: The square of its 
Minkowski norm. 

The electromagnetic field strength corresponds to the bivector F = E + iB, 
and Maxwell's inhomogeneous equation (11.21) translates to 

VF = -A-kJ. (2.4) 

By operating with V on both sides, this equation leads to 

VVF = -4vrV J, (2.5) 

where VV is a scalar operator: The d'Alembertian. Therefore, since F is a 
vector, and because operating with a scalar operator on a vector can only lead to 
a vector, taking the scalar part on both sides of (12.51) automatically leads to (12.31) . 
the biquaternion form of the tensor identity (11.11) . Consequently, in biquaternion 
formalism the necessity of the local conservation of J is even more transparent than 
in tensor formalism: It is immediately seen that this conservation is a consequence 
of the formal structure of electrodynamics, and for this reason more fundamental 
than, for example, gauge invariance flU p. 676]. 

Finally, to end this brief introduction to quaternion methods, we give the 
definition of F in terms of the four-potential A = ip — iA, i.e., 

F:=VaA, (2.6) 

where the the binary operator A as in Q A R means taking the vector part of the 
product QR. 
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2.1 Paul Weiss's causal coordinate system 

The most important covariant null-quantity of field theory is, in tensor notation, 
the four- vector interval = — between the position Z^ of moving particle 
and the point at which the fields associated with that particle are observed. 
Indeed, the nullity of the Minkowskian distance 

RpB? = \x - z\ 2 - {t x - t z f = 0, (2.7) 

means that the points X and Z are positioned in such a way that a field observed 
at X can be causally related to its source at Z, and vice- versa that a particle at Z 
can be causally influenced by the field generated by another particle at X. In the 
first case causality is implied by the condition that the time t z is preceding t x , i.e., 

t z < t x . (2.8) 

The corresponding root t r := t z of equation (12.71) is customarily called the 'retarded 
time,' although 'causal time' would be a better word to suggest that it is the time 
t r at which the source located at z produced the effect that is later observed at the 
location x at the time t x . 

In biquaternion notation the observer's four- vector is X = it x + x and that of 
the source Z = it z + z. Equation (12.71) is then written 

RR = (X - Z) (X - Z) = 0. (2.9) 

If the world-line of the source is expressed in parametric form as a function Z(r) 
of the proper time r = jt z , where 7 = — ft 2 is the Lorentz factor and 

(3 = v/c the three-dimensional velocity, its four- velocity Z is given by 

Z(t) =7(1 -0), so that ZZ=1, (2.10) 

where the dot corresponds to differentiation with respect to imaginary proper time, 
i.e., Z = dZ/dir. 

Since biquaternions form an algebra, the four- velocity has a unique decompo- 
sition as a product 

Z(t) = BB + , suchthat EB=\. (2.11) 

Here the operator () + denotes quaternion biconjugation, namely the operation 
(s + v) + — (s + v) = s* - v* such that (QR) + = R + Q + , which combines 
quaternion and complex conjugations. Eq. (12.1 II) is the spinor decomposition of 
the four-velocity, and the unit biquaternion £>(r) biunivocally defines (up to a 
factor ±1) the general Lorentz-boost to the frame in which the arbitrarily moving 
particle at the point Z{f) is instantaneously at rest. 
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In biquaternions the causal relationship (12.91 ) between the points X and Z 
can be expressed by means of the spinor decomposition of the interval X — Z, 
discovered in 1941 by Paul Weiss ifTOll . i.e., 

X - Z = R = £B(i + v)B + . (2.12) 

Here £ is an invariant scalar and V = i>(9, 0) a unit vector so that i + u, and 
consequently R, are null biquaternions. Equation (12.121) is therefore an explicit 
parametrization of R in terms of four variables: The invariants £ and r, and the 
two angles 9 and <p characterizing the unit vector V in the instantaneous rest frame 
of the particle. Thus, in that frame, i.e., when B = 1, the interval X — Z reduces to 
the quaternion £(i + z7), which shows that in that frame the vector £z7 corresponds 
to the ordinary radius vector x — z, and that the distance £ appears to an observer 
at rest in that frame as the ordinary distance \x — z\. On the other hand, when 
B 7^ 1, Eq. (12.121) provides a general parametrization of the null-interval R, that 
is, in geometrical language, of the 'light-' or null-cone originating from Z. 

In the case where Z, £, u(9, 0), and B are evaluated at the retarded proper time 
r r := 7t r , Eq. (|2.12l) defines a causal coordinate system erected at the retarded 
four-position Z(r r ), i.e., 



X(T,£,9,(f>) = Z + ^B(i + P)B- ] 



(2.13) 



Paul Weiss called it 'retarded coordinate system,' but we prefer the term causal 
coordinate systemQ From the form of equations (12.121) and (12.131) it is immediately 
seen that £ can be expressed in terms of R and Z, namely as the scalar product 



ii = ZoR 



(2.14) 



which is the usual definition of the retarded distance £. As £ is equal to the 
difference ^{t x — t z ) = r — r r , we also have r r = r — £. 

In the present paper the spinor representation (12.121) of R will not be used 
anymore. Instead, to make all £ dependencies explicit, we introduce the 'unit' 
null-quaternion K such that 



K(0,4>):=R/£ 



(2.15) 



where as above the condition r = r r recalls that all quantities are evaluated at the 
retarded proper time r r . In the following, for simplicity, this condition will be 
specified only for the main equations. 



'Equations (12.12b and (12.131 1 describe the motion of a generally accelerated particle from the 
standpoint of an inertial frame that is momentarily at rest with respect to the particle. In such 
a case the actual comoving frame of reference of the accelerating particle is noninertial, which 
implies that the object of anholonomity has to be introduced to properly define differentiation ifTTI . 
In this paper this is taken care of by an intrinsic definition of the causal derivatives introduced in 
subsection (2.3) below. 
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2.2 Kinematical invariants 

In this paper we shall have to calculate expressions containing up to the fourth 
proper time derivative of Z. This implies that apart from £ there will be three more 
invariants, namely the acceleration, biacceleration, and triacceleration invariants 
defined as: 

k:=ZoK, (2.16) 
i x ■- ZoK, (2.17) 

(:=ZoK. (2.18) 

In the literature these invariants are sometimes defined as Z o R, Z ° R, etc. As 
R = £K this multiplies our invariants by a factor £, which is better left explicit 
since £ is a coordinate (even though the combinations re£, x£ and arise very 
often during calculations because of dimensionality). 

We shall also use a number of well-known kinematical identities that are easily 
found by iteratively derivating Eq. (|2.10l) . and which yield the invariants: 



ZoZ = ZZ = 1, (2.19) 

ZoZ = 0, (2.20) 

Zo~Z = —ZZ = -A 2 , (2.21) 

ZoZ = -3Z o Z, (2.22) 



where the invariant A is the magnitude of the acceleration. 

2.3 Causal derivatives 

When differentiating quantities depending on X and Z with respect to the coor- 
dinates of the point of observation X, one has to include the partial derivatives 
relative to both X and Z since they are related by the condition RR = 0, which 
insures causality. Because Z = Z(r r ) is an implicit function of r this can be done 
with the help of a simple rule obtained by differentiating Eq. (12.91) . i.e., 

2(dX -dZ)oR = 2(dX - dir r Z) oR = 0. (2.23) 

Thus, 

dir r = dXo -JL= = -idX o K, (2.24) 
ZoR 
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where we used (12.141) and the definition of K. But, by definition of the total 
derivative, 

dir r = dX o Vzr r , (2.25) 

where V is the four- gradient operator (|2.2I) . Therefore, comparing (12.241) and 
(l2~25l we find 

Vir r = -iK. (2.26) 

Consequently, for an expression E = E(X, r), where the argument X corresponds 
to an explicit dependence on X, and r to the proper time, we have the covariant 
differentiation rule 



VE(X, t) = VE(X) - iKE(r) 



(2.27) 



We now use this rule to derive a number of simple relations which will be 
useful for calculating more complicated expressions in the following. The first 
one is 

VR = 4 + iKZ. (2.28) 

where we used VX = 4 which is obvious. Further we have 

N := V£ = (1 — k£)K - iZ, (2.29) 

which defines the normal N to the proper tube of radius £ surrounding the world- 
line, and which has the properties 

jVojV=l-2<, (2.30) 

NoZ = -iK£, (2.31) 

NoZ = k(1-k£), (2.32) 

N o'z = iA 2 + -Ox- (2.33) 

Using the same method, i.e., the differentiation rule (12.271) and the previously 
obtained results, we also find 

VK = -VVr r = -, (2.34) 

ViV = VV£= (1-2k0|. (2.35) 
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Similarly, we have the four-gradients 



Vk£ = Z + x^K, (2.36) 
V X £ = -%'Z - {A 2 + tf)K. (2.37) 

The second of these involves a fourth proper-time derivative because C, = Z o K . 
Accordingly, if we calculate the d' Alembertian of k£, i.e., the gradient of equation 
(|2.36l) . we expect such a fourth derivative to arise. However, 

VV/cf = 4*. (2.38) 

The reason for this is that when calculating the gradient of the second term in 
(12.361) we get 

V X £K = {-it - {A 2 + ()K)K + xZVK, (2.39) 

so that the term containing ( disappears because KK = 0. As will be seen, this 
is a general feature: Causality insures that no derivatives of Z beyond the third 
appear in any physically meaningful quantity considered in this paper. 



3 The customary formulation 

The potential of the Lienard-Wiechert field can be obtained by a number of methods 
that are explained in many text books, e.g., |[T2l[T3l[l4l[T5l[T6l. Most frequently 
it is obtained by working in the Lorenz gauge, and by means of a Green's func- 
tion assuming that the point-charge can be represented by a three-dimensional 
5-function. In biquaternion notation this source current-density can be written in 
the following covariant form 

j s = ^ao, (3-D 

whose normalization corresponds to the global (or integral) form of charge con- 
servation, i.e., 

r-l'K pit poo 

j d(j) d6 sin6 / dH 2 J = eZ. (3.2) 
Jo Jo Jo 

The resulting potential has the remarkable simple form 

(3.3) 

T=T r 
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from which, one can calculate the electromagnetic field strength bivector F = 
E + iB according to its definition, Eq. (12.61) . which can be rewritten in equivalent 
form as 

F = V AA = VA-VoA = B -ip, (3.4) 

where 

B := VA, and the scalar ip : = V o A, (3.5) 

is the invariant which is set equal to zero in the Lorenz gauge. 

Applying the rule (12.271) and using equation (12.291) it is easily found that 

of which the scalar part is, using equations (12.141) to (12.161) and being careful not 
to subtract potentially infinite quantities, 

^ LW = -e((l - - ^ 1 ^ 1 ) T=Tr = 0, 0, (3.7) 

which shows that the Lorenz gauge, i.e., ip = 0, is only satisfied when £ ^ 0. 
However, since the electromagnetic field strength is a vector, the fact that the scalar 
ip LW is undefined at £ = has no direct consequence on the field, even though 
the field F LW is derived from the same potential A LW , which just like F LW is 
undefined at £ = 0. Thus, taking the vector part of (13.61) the Lienard-Wiechert 
field strength is 



= -e[(l - «{)_ -i—^j^. (3.8) 

Using again the rule (12.271) we can calculate the corresponding charge current 
density according to (12.41) . However, since ip LW = for £ ^ 0, it is easier to use 
(|3~4|) and rewrite (f2~4|) as 

jLW = _ V vA LW , V£ ^ 0, (3.9) 

which, in fact, is the standard form used in the customary formulation assuming the 
Lorenz gauge. Then using for VA LW equation (13.61) . we get after some calculations 

^ LW = tMj. ~ ~ a " 2 ^)(|r " <3-i0) 
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In deriving this result we have been careful not to subtract potentially infinite 
quantities in order to see how (^-functions, such as in the current density (13.11) . 
arise. Indeed, when £ 7^ the current (13.101) is zero, but since it is undefined for 
£ = one has to substitute an appropriate 5-function for the ill-defined differences 
— l/£ n . The standard prescription for doing this is to use global charge 
conservation, i.e., equation (13.21) . which implies that the expression 2/£ 3 — 2/£ 3 
should be replaced by — 5(£)/£ 2 - But there is no such constraint for the l/£ 2 — l/£ 2 
expression in the first term of (13.101) . However, if its origin is traced to the l/£ 
factor in the potential, one finds that it should be to be consistent with 

the 2/£ 3 — 2/£ 3 expression. With this prescription we get 

(3.11) 

T = T r 

This current density differs from (13.11) by the presence of two additional terms 
which depends on the acceleration. However, when integrated over the whole 
three-space as in (13.21) . J LW yields the same total current eZ as J s , because after 
multiplication by the volume element d<p dO sin 8 £ 2 the acceleration dependent 
terms does not contribute to the radial integral since £5(£) = 0. Unfortunately, if 
one tries to verify that the charge-current density J LW is locally conserved, and 
consistently uses again the differentiation rule (12.271) . one finds 

This result is non-zero, and distributionally well-defined in the sense that the 
measure J d£ £ 2 <9 M J^ w is finite. The current density J LW is therefore not locally 
conserved, because according to (11.11) the divergence (13.121) should be identically 
zero. 

Moreover, if instead of J LW one tries to verify the local conservation of the 
customary source current-density J s defined by (13.11) . one also finds that it is not 
locally conserved. Indeed, as this current is very simple, this is immediately seen 
by calculating its divergence, which consists of three terms, i.e., using equations 
tiTZh . (I239T) . (I2~T61) . and (OH) 

^VoJ s = -*Zo(z)i5(0 
+ FoZ(— )*(£) 

= (-iK + 2ik + ^ 0. (3.13) 



J 



LW 
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Only in the inertial limit, that is for non- accelerated motion, are the current 
densities J LW and J s locally conserved. This means that something is wrong 
in the customary formulation of the electrodynamics of an arbitrarily moving 
point-charge, or else that something is inconsistent and needs to be clarified. 
For instance, to avoid any misunderstanding, it should be recalled that there is 
large class of conditions under which a simple current density like (|3.1I) is locally 
conserved. These comprise the case where instead of being related by a causality 
condition such as (12.91) . the X and Z variables are independent, e.g., lfT3l p. 72] 
and EH p. 139]. 

4 The locally conserved charge-current density 

In this section we derive the locally conserved current of an arbitrarily moving 
point-charge by making no other assumption than that its potential reduces to a 
scalar function cp(£) in its instantaneous rest frame. Thus, instead of taking for 
the potential the Lienard-Wiechert form, i.e., equation (|3.3I) . we will assume the 
general form 



and seek under which conditions it leads to a conserved point-like current density. 

This method is motivated by the observation that any current distribution of 
the general form 



where S and T are any smooth four-vector functions, will satisfy global charge 
conservation because of the distributional identity = 0. We therefore 

suppose that the absence of local current conservation could be due to an incorrect 
handling of the singularity at £ = 0, which should in fact lead to a current having 
a more complicated form than equations (13.11) or (13.1 II) . Thus, to proceed step 
by step and be fully general, we replace the 1 /£ factor in the potential (13.31) by a 
function (p(£) in 04. 1 1) that is finite and indefinitely differentiable except possibly at 
£ = 0. This will lead us to a differential equation for a conserved current, of which 
we can study the solutions in terms of both regular functions and distributions. 

Moreover, in order to make gauge invariance explicit, we write Maxwell's 
equation for the current density J directly in terms of A, i.e., from (|2.4I) and (13.41) . 



A := eZ(r)<p(0 



(4.1) 




(4.2) 



VVA - VV> = -471- J, 



(4.3) 



where ip is the invariant already defined in (13.51) . 
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4.1 General field-strength and gauge invariance 

To find the field strength (|3.4I) we first calculate B, which using (12.291) and (12.271) . 

is 

B = VA = VeZip = eNZ<p' - eiKZ<p. (4.4) 
Using (12.311) and (12.161) . its scalar part is 

ip = —ein^ipi, (4.5) 

where 

Pi ■■= <f' + 7<f- (4-6) 

The field strength, equation (13.41) . is then obtained by subtracting the scalar ip from 
the biquaternion B = V A what is equivalent to taking its vector part, i.e., 

F = V A A = e~N A Zip' - eiK A Zip, (4.7) 

which is gauge invariant. Indeed, by adding to A the four-gradient of any scalar 
function f(X, Z) of the spatial coordinates and time, i.e., by making the substi- 
tution A -> A + V/, we see that B -> B + VV/ while ^ -> V + VV/, so 
that F — i> F and consequently neither F nor J depend on the choice of the gauge 
function /. 

Moreover, as ip has been subtracted from B = VA, the field strength F does 
not depend on ip nor on cpi, only on ip and cp' separately, as can be seen in (14.71) . 

4.2 General charge-current density 

To find the general current density (14.31) . we begin by calculating VB, what is 
simpler to do if it is rewritten as eVV 'Zip. Thus, 

WB = e ( V VZj <p + e (v<p) (vz) + eV (V<p) Z + e (yV<p) Z, (4.8) 

where we have used the 'mobility' of the scalar function ip to put it at the right 
place for the differentiations which are made within the ranges indicated by the 
parentheses. The four terms in this expression are 

(VVZ) <p = -i (vKZ^j v = -2i-Z<p, (4.9) 

(Vp) (VZ\ = (Nip')(-iKZ) = -ZKZip', (4.10) 
v(V<pjz = -iK(N(p')Z = -K~Z~Zip', (4.11) 
(VV^JZ = (ip'VN + NN<p"^Z = (1 - 2Kg)((f/' + -(f/)Z, (4.12) 
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where in the course of the calculations the identities of section [2] were used, 
including the relation KK = which enabled several simplifications. Combining 

these four terms, and noting that ZK + KZ = 2Z o K = 2i, we find 

Vfi = -e2iZ Vl + e(l - 2n^)Z(p 2 . (4.13) 

Here, 

<P2(0 = <f/' + l<S, ( 4 - 14 ) 

has been introduced as an expression distinct from (pi, even though it could be 
expressed in terms of cpi as cp 2 = f[ + <fi/£- This is because we want to avoid 
any subtraction of potentially undefined terms, and therefore keep everything 
expressed in function of ip and its derivatives as they arise. Moreover, as seen 
in (14.121) . (f2 arises in the expression of the d'Alembertian of the scalar function 
which corresponds to the field in spin-zero electrodynamics, and is therefore 
of interest on its own. 

We next calculate V^>, which according to (14.51) is 

V-0 = -ei^iV/e£ - etK£N(f/ v (4.15) 

Using (12.291) and (12.361) this can be rewritten as 

W = -ei(Z + x£K>i - ei<((l - k£)K - iZ)^, (4.16) 

which shows that whereas VB given by (|4.13l) is function of Z and Z, Vip, and 
thus the current density J given by (14.31) . are moreover function of Z through the 
invariant \- 

Finally, combining (14.131) and (14.161) . the general form of the charge-current 
density is 

-4vrJ = + e(l - 2k£)Z<P2 
- ei(Z - x£K)ip x 

+ ei((l- k£)K -iZ)K£<p[. (4.17) 

Comparing with (|4.12l) . the first line is seen to come from the d'Alembertian of 
(p(£), which is also at the origin of ^(O- On the other hand, the second and third 
lines, as well as <£>i(£), come from both V B and Vip. 

In the inertial limit, Z = Cst, this current-density reduces to the simple 
expression 

-4vrJ = +eZip 2 , (4.18) 

which when </? = l/£ is just the Green's function usual source current density 
(13.11) . On the other hand, the Lienard-Wiechert current-density (|3.1 II) is obtained 
by ignoring the contributions coming from ip, that is by taking equation (14.131) 
equal to —An J and assuming cp = l/£. 
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4.3 Divergence of general charge-current density 



Since we have the general form of the charge-current density we can now explicitly 
verify that V o J = 0, i.e., equation (11.11) . Of course, we expect that to be true if 
there is no algebraic error in our calculations. But it is nevertheless important to 
do it since we want to make sure that nothing special arises in terms of handling 
ill-defined quantities, and also because we want to know under which precise 
conditions the charge-current density is locally conserved. 

However, since the calculation of V o J is tedious, and basically consists of 
repetitively using the previously used methods and identities, we will just give the 
result, which is 

V o J = + 2eix^(i(f[ + Tfi) - V2) 

+ em£(l - 2<) (V; + - (<p' 2 + ^2)) . (4.19) 

As can be seen, this expression contains only the three invariants £, k, and 
X- There is therefore no dependence on the fourth proper-time derivative of 
Z(t), despite that the initial expression, the potential (14.11) . depends on the first 
derivative through Z and z£ = R o Z. As seen at the end of section [2l this is 
because causality insures that no derivatives of Z beyond the third appear in any 
physically meaningful quantity considered in this paper. 

On the other hand, concerning <£>(£), the condition of validity of equation 
(14.191) is, as expected, that it is three times differentiable, provided £ 7^ 0. Hence, 
if we assume such a function for y?(£)> an d calculate the expressions in the two 
big parentheses, we find that they are both identically zero, implying that the 
charge-current density is indeed conserved for any £ except possibly for £ = 0. 
However, if is a distribution rather than a regular function, which means that 
the derivatives of y?(£) will be defined even at £ = 0, the current density J will 
also be locally conserved, but the condition £ 7^ will no more be required. 

Finally, in order to illustrate the importance of not having made any simplifi- 
cation before getting V o J in its final form as in (14.191) , we remark that the first 
line of this equation would not be there had we expressed ip 2 in terms of ipi, a 
possibility that was mentioned below equation (14.141) . 



4.4 Locally-conserved general charge-current density 

We finally come to the Lienard-Wiechert case, and therefore specialize to <p = l/£. 
For ^ we have then 

^i(O = -^ + ^=0, and MO =+^-^=°> ( 4 - 2 °) 
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which imply that all terms in (|4.17l) are zero, including £ip[ = £cp 2 — (fi. The 
current density J is then everywhere zero, except at £ = where it is undefined. 
We therefore interpret y?i(£) as a distribution, and use the theorem stating that a 
distribution which has its support only in one point, say the origin, is a linear 
combination of the 5-function and its derivatives up to a certain order [2, p. 784], 
El p. 443]. Thus 

V£ > 0, </?i(0 = -M), (4-21) 

which because of dimensionality comprises a single 5-function, and whose nor- 
malization will turn out to be consistent with (13.21) . Similarly, we also interpret 
<y9 2 (£) as a distribution, and, since by (14.141) it can be expressed in terms of <^i(£)> 
we ge|§ 

V^O, ^) = ^ 1 + I^ 1 = -1^). (4.22) 
It remains to substitute (14.211) and (14.221) into (14.171) . and the locally-conserved 



(4.23) 



charge-current density (14.31) is finally found to be 



e (Z .Z + 2kK ,. n „ wAr/.x 



e /Z .Z + 2kK 
This result leads to several observations: 



1. The current density J is much more complicated than the simple current 
(|3.1I) : It depends directly on the three invariants £, k, and x, as wei l as on 
the two four- vectors Z and Z; indirectly on the biacceleration Z through the 
invariant x', and, finally, on the angular variables through the null four- vector 

K(9,<P). 

2. The dependence on the third derivative of Z is consistent with the Lorentz- 
Dirac equation and with the Schott expression of the self-force, in which Z 
also appears, because the self-interaction force involves a product of J with 
the self-field. 

3. Equation (14.231) has the most general distributional form of (|4.2I) . in accord 
with the theorem cited above equation (14.211) . 

4. With (pi given by (14.211) the invariant ip defined by (14.51) is 

i/j = -eiK5(£). (4.24) 



incidentally, comparing with the discussion above equation (13.111 1. we have rigorously justified 
the prescription that led to the non-locally-conserved charge-current density (13.11b . 
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Thus, the gauge can be the Lorenz gauge ip = only for £ ^ when 
i.e., when the acceleration is non-zero. 



5. The equation y?i(£) = y?' + y?/£ = 0, V£ ^ 0, has only one solution: 
This singles out the corresponding potential as being the only one such that 
the current density of a point-charge is conserved and thus given by (14.231) . 

5 Straightforward derivation 

While the previous section's derivation of the locally-conserved current of an 
arbitrarily moving point-charge is rigorous, it has one important shortcoming: The 
derivation basically consisted of using distribution theory to find the correct form of 
that current, without having explicitly specified the forms of the potential and field 
leading to it, even though we know from section [3] that something is deficient in 
their customary formulation since they do not lead to the locally conserved current. 
Indeed, what we expect from a consistent application of Maxwell's theory is that 
once we have properly formulated the potential or the field, the current should 
derive from either of them in a straightforward manner, that is without having to 
'mend' some ill-defined result to get the correct one. 

From a mathematical perspective, this shortcoming consists of having con- 
sidered just (fx and (p 2 as distributions, whereas according to equations (14.61) and 
(14.141) they both derive from ip which appears in the potential (14.11) already, and 
which therefore should also be defined in the context of distribution theory. In 
other words, we have not looked at the origin of the ^-functions arising in cp\ and 
(p 2 , which according to Schwartz's structure theorem of distribution theory must 
arise from the partial differentiation of some continuous function. In fact, the 
generating function leading to these (^-functions is easily found. This is because, 
in three-dimensional notation, the retarded distance Eq. (12.141) reads 



where p is the unit vector in the direction of x — z. The retarded distance is 
therefore proportional to an absolute value, and for that reason has a discontinuous 
derivative when x — > z, i.e., at £ = 0. Thus, when calculating derivatives of func- 
tions expressed in the coordinate system {£, 9, 0}, one must carefully distinguish 
between the coordinate £ and the absolute value |£|, a crucial observation that was 
first made by Tangherlini [|7] p. 511-513]. 

Consequently, as is explained in details in references [ffl[5l|6]|, the potential of 
an arbitrarily moving accelerated point-charge must be written 



£ = |z-z|7(l 



(5.1) 



^4 = 4 T (0 



(5.2) 



18 



where T(£) is the generalized function defined as^| 

{undefined r < 0, 
r = 0, and -f T(r) = <5(r). (5.3) 

ar 
+1 r > 0, 

When the definition (|2.6I) and the causal differentiation rule (|2.27l) are now used 
to calculate the field strength starting from the potential (15.21) . the corresponding 
current density (|2.4I) is directly found to be the conserved one, i.e., equation (14.231) . 
However, instead of the customary Lienard-Wiechert field, equation (13.81) . the field 
strength is now 

(5.4) 
(5.5) 

which apart from the presence of the T-function multiplying F LW on the first 
line, has an additional 5-like contribution on the second. Indeed, starting from the 
general form of the field F given by (14.71) . and taking for (p the function implied 
by the potential (15.21) . it suffice to substitute 

¥>=~T(0, and ^' = -I T (0 + ~5(0, (5.6) 

in (gJl) to obtain (T54M531) after use of the identity N A Z = (1 - A Z. 

Since both the T-factor and the 5-like contribution in F are necessary to obtain 
the current density satisfying the local conservation identity (11.11) . it is clear that 
the customary F LW cannot lead to such a current. In fact, it is by calculating the 
current density immediately from the potential as in (14.31) — that is by ignoring 
that the field strength could be different from the customary one — that after many 
unsuccessful attempts the author discovered the conserved current density in July 
2003. 

Indeed, the <5-like field on the second line of (15.4115.51) is absolutely necessary 
since it carries an essential part of the information about the nature of the singularity 
at £ = 0, which is not simply due to the divergence of l/£ when £ — > 0, as in the 
customary formulation, but to that of 1/|£| where the absolute value leads to the 
expression T(£)/£ in the potential (15.21) . 

As this T-function in the potential (15.21) is at the origin of the ^-functions 
arising in the conserved current, we can come back on the functions (fi and ip 2 

3 Intuitively, T(r) can be seen as equivalent to the sign function sgn(r) for r > 0. 
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+ e((l-<) 



K A Z _ .K A Z) 
K A Z 



T(£) 



defined in the previous section, and observe that starting from = T/£ we have 
(fi — 5/i and (f 2 = 5'/£. Comparing to equations (|4.21l) and (|4.22l) . we also have 
ifii = —5' D and (f 2 = —5'^, where 5r> is the usual 'Dirac 5-function' which has the 
property 5' D = — <5d/£ characteristic of Schwartz's distribution-theory. However, 
since ip\ and if2, and thus ip, enter the conservation equation (14.191) — which is an 
identity — that property is not necessary for insuring local charge conservation. 
This means that any generalized function such that T' = 5c (where 5c ^ #d is for 
example a Colombeau generalized function [6J) is sufficient to have local charge 
conservation provided 5' c and 5 C exist. 

The major implication of the necessity of the ^-contribution (15.51) to get an 
identically conserved charge-current density is that the field (|5 .4145 .51) cannot be 
interpreted as a distribution in which the 5-contribution (15.51) can be ignored if the 
field is twice differentiated to calculate the divergence (11.11) . The reason is that 
every successive differentiation introduces new terms which must be retained in 
order to get a result that is identically zero The field (I54U531) and the current 
density (14.231) must therefore be interpreted as nonlinear generalized function, such 
as Colombeau functions [6], in which all terms can be or can become significative 
in further calculations, even though none of the nonlinear properties of such 
functions have to used in the present paper. 



6 Discussion 

In this paper we have derived the proper formulation of the potential of an arbitrarily 
moving point-charge, equation (|5.2I) . which leads to the current-density (14.231) that 
is locally-conserved, and to the field strength (15.41) which contains an T-function 
factor and a 5-function term that are absent in the customary Lienard-Wiechert 
form of that field. Since the reasoning leading to these results is fully general, and 
therefore valid in the inertial limit, the Coulomb potential in the rest frame of a 
point-charge should thus be written ||4l[5l[6]| 

<p c (r) :=e-T(r), (6.1) 
r 

which implies that the corresponding Coulomb field 

E c (r) = -Vfc = e^T(r) - e^5(r), (6.2) 

also has a 5-function term that is absent in the standard formulation of that field. 
As shown in [5J, this <5-term implies that the self-energy of a point charge is not a 

4 This can be illustrated by the following example: the function /(£) = £ _1 <5(£) yields zero 
when evaluated on any test function T(£) in R 3 . On the other hand, /'(£) = — 2£ -2 £(£) yields 
— 2T(0) / so that differentiation and evaluation on test functions do not commute in general. 
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1/r distribution but a 5 2 {r) singularity, so that the self-energy of a point charge is 
fully concentrated at r = 0. 

As the necessity to supplement the customary Lienard-Wiechert and Coulomb 
fields by a 5-function term is quite unexpected, it is important to stress that this 
conclusion was reached while remaining entirely within the realm of standard 
Maxwell theory, that is without making any modification to that theory. The 
only thing that was done was to carefully analyze the nature of the singularity at 
the position of the point-charge, and to use distribution theory to formulate the 
implications stemming from it. Moreover, while the derivation of the conserved 
current was somewhat indirect in sectional its straightforward derivation from the 
potential 



in section \5\ and the rigorous mathematical justification of that potential, make 
that the replacement of the Coulomb potential and field by the expressions (16.11) 
and (16.21) is of absolute physical necessity. 

Another aspect of the present paper may look somewhat surprising: distribution 
theory has since long been applied to classical electrodynamics by a number of 
researchers, and no modification to the Lienard-Wiechert potential and field was 
found to be necessary, see, e.g., lfT71[T8l . The reason is that these researchers did 
not question the standard form of the Lienard-Wiechert potential and of the fields 
deriving from it: they simply postulated that they were distributions in order to 
consistently work with their singularities at £ = 0. On the other hand, what we 
have done in this paper was to uncover the proper form of the potential leading to 
a point-charge current-density that is locally conserved — which implied that just 
like in nonlinear generalized function theory all terms obtained by differentiation 
had to be kept until the end of the calculation in order to get an identity. 

In doing so we have in fact discovered more than we had hoped for, namely a 
striking result that can be phrased in the form of a unicity theorem: 

Theorem [Unicity of the four-potential of an arbitrarily moving 
point-charge] The only potential of the form 



^4 = 4 T (0 



(6.3) 



A = e<p(OZ(r) 



(6.4) 



where (p(£) is a scalar function, such that the charge-current-density 



Ait J = V(V A A), 



(6.5) 
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is a locally conserved 5-like nonlinear generalized function, i.e., such 
that 



V o J = 0, Vf , (6.6) 
AO =0, V£ ^ 0, (6.7) 



is the unique solution of the equation 



namely 



<p'(0 + 7*>(0 = £*(0. (6-8) 



AO = 7 T (0, (6-9) 



provided all terms obtained by differentiation are kept until the end of 
the calculation. 



Finally, since the formulation of the potential and field of an arbitrarily moving 
point-charge presented in the present paper is new, there is an apparent contradic- 
tion with the fact that the customary Lienard-Wiechert formulation is an agreement 
with so many applications of classical electrodynamics. 

There is however no contradiction, since, on the contrary, the results of this 
paper are in full agreement with the fundamental principles of electrodynamics 
and mechanics. For instance, if the conserved current (14.231) is introduced in a 
Lagrange function as the scalar product J o A ext with the potential of an external 
field A cxt E C°°(M 4 ), the differences between that current and the current J s 
defined by equation (|3.1I) have in general no influence since they disappear upon 
integration over the whole space. The same is true for the derivation (which also 
involves an integration) of the potential of an arbitrarily moving point charge, i.e., 
equations (13.31) or (15.21) . by means of a Green's function. 

Thus, the principles of physics imply that the position Z and velocity Z of a 
point-charge are sufficient to determine the potential of its field, while the precise 
formulation of that potential given by (|5.2I) is necessary to determine the complete 
field and conserved current-density, which include terms that are function of Z 
and Z. In other words, while J s is sufficient as a source to determine uniquely the 
potential of an arbitrarily moving point-charge, the conserved current J deriving 
from this potential can be very different from J s . Indeed, without the factor e, J s 
is merely a velocity distribution associated to the world-line of the point-charge. 

In conclusion, the formulation presented in this letter will make little difference 
in most engineering-type applications of classical electrodynamics. Indeed, as can 
be seen by studying a number of examples, the instances in which the full details of 
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the current density (14.231) are strictly necessary, and the additional ^-contributions 
in the fields (15.4145.51) and (16.21) are essential, include fundamental problems like 
calculating the interaction of a point-charge with itself [0, and similar nonlinear 
problems in which classical electrodynamics is apparently not consistent. For 
example, the self-force on an arbitrarily moving point-charge is given by F^ U J U 
where F^ v is given by (15.4145.51) and J u is the full conserved current density (14.231) . 
i.e., such that d u F fMU = —AnJ^, not the customary charge-current density (13.11) : a 
frequent confusion which is a cause of known difficulties. The resolution of such 
internal contradictions is the subject of several forthcoming publications lfT9ll20l . 

Finally, we stress that the methods used in present paper to derive the potential, 
field, and charge-current density of a point-charge can be extended to the derivation 
of those of an electron, i.e., a point-charge endowed with a magnetic dipole moment 
ft and spin 0, whose four-potential expressed in the spinor notation of section [2] 

is eh 

^electron = B(~ - Z^)^ + T(0. (6.10) 
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